Radiative falloff in the background of rotating black holes 
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We study numerically the late-time tails of linearized fields with any spin s in the background of a 
spinning black hole. Our code is based on the ingoing Kerr coordinates, which allow us to penetrate 
through the event horizon. The late time tails are dominated by the mode with the least multipole 
moment I which is consistent with the equatorial symmetry of the initial data and is equal to or 
greater than the least radiative mode with s and the azimuthal number m. 
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The late time dynamics of black hole perturbations 
has been studied for over three decades. Complete un- 
derstanding of the late-time dynamics is available for a 
Schwarzschild background: Generic perturbation fields of 
either scalar, electromagnetic, or gravitational fields de- 
cay at late times along an r = const curve as an inverse 
power of time. Specifically, linearized fields (the scalar 
field itself, or the Teukolsky function in the gravita- 
tional case) decay as (assuming that the initial 
data have compact support and are not time-symmetric), 
where £ is the multipole moment of the perturbation field 
Q, 1^ 1^ . This behavior was confirmed also for fully non- 
linear collapse of spherical scalar fields 0, 0| . The mech- 
anism which is responsible for this behavior is the scat- 
tering of the field off the curvature of spacetime asymp- 
totically far from the black hole. 

Because it is only the asymptotically far geometry 
which determines the behavior of the late-time tails, it 
is natural to expect similar behavior also when the black 
hole is rotating 6]. Because spacetime is not spherically 
symmetric, however, spherical-harmonic modes do not 
evolve independently. Specifically, taking the initial data 
of the perturbation field to be a pure mode, other 
modes are excited. Intuitively, all the modes which are 
not disallowed [by symmetry requirements (such as the 
equatorial symmetry of the initial data) or dynamical 
considerations (such as that only modes with —£< m < £ 
are allowed)] will be excited. In particular, modes with 
£ values smaller than the original £ will be excited, and 
will dominate at late times. (Notice that because the 
background is axially symmetric, modes with different 
values of m are not excited when linearized perturbation 
theory is applied.) Accordingly, the late-time dynam- 
ics is dominated by the mode with the least £ which is 
excited, namely the smallest £ which is not disallowed. 
That is, all modes £ which are not smaller than |r7i| and 
|s|, where s is the spin weight of the field, and which re- 
spect the equatorial symmetry of the initial data, will be 
excited. The falloff rate is then ^-"(2^min+3)^ where i^min is 
the smallest mode which can be excited. 



Despite the simplicity of this intuitive picture, recent 
papers report conflicting results. An analytical analysis 
by Hod — in which the author attempted to find the 
asymptotic behavior of the fields in the spacetime of a 
Kerr black hole — yielded results which are more com- 
plicated: The decay rate for a scalar field is predicted 
by Hod to be t-(2^*+3) H £* = m or £* = m + 1, 
^-(r+m+i) if i* - m > 2 is even, and t-(^*+™+2) if 
£* — m > 2 is odd, where £* is the initial value of 
£. For gravitational perturbations Hod's formula is 
f-ie +io+3-q) j-j-Qj- axisymmetric perturbations), where 
£q is the radiative mode with the least value of £, and 
q — min(^* — £o, 2). [Different, apparently conflicting re- 
sults were reported by Barack and Ori 9] . Those authors 
assumed that the mode £,m = is present in the initial 
data (for s = 0), as a result of which it is not straight- 
forward to confront their predictions with Hod's.] 

Although Hod's results could be relevant for an in- 
termediate regime for carefully chosen parameters, they 
make only little sense for describing the intended asymp- 
totic late-time behavior. These eerie conclusions imply 
that some sort of a "memory effect" takes place: the 
field somehow "remembers" its initial configuration, de- 
spite being a linearized field. We do not believe that such 
a "memory effect" is reasonable: Take the initial data at 
the time to to be those of the pure mode £* , such that £* 
is significantly larger than ^min- At the time ti > to the 
field includes, in addition to the mode £*, also contribu- 
tions from modes £ < £* because of the excitation of other 
£ modes. Now the fields a.t t — ti can be construed as the 
initial data of a new evolutionary problem. In the new 
problem the initial data are a mixture of modes, such 
that modes £ smaller than £* are present 0. Because 
the mode with the smallest existing £ value dominates 
at late times and determines the decay rate of the tail, 
we can see no way in which the £* mode can determine 
the asymptotic late-time tail, unless £* determines which 
modes can and which modes cannot be excited. As in 
the spacetime of a Kerr black hole it is hard to see how 
a scenario in which modes which arc not disallowed can 
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still be excluded, we conclude that "memory effects" are 
not to be expected. Hod's results, if correct, suggest to 
us that an hitherto unsuspected mechanism of selection 
rules inhibits the excitation of otherwise allowed modes. 
Such a counter-intuitive theoretical reasoning must have 
strong numerical support in order not to be discarded. 

Conclusions which apparently are similar to Hod's were 
obtained more recently by Poisson j^] , who analyzed the 
scalar-field tails in a general weakly-curved, stationary, 
asymptotically flat spacetime. We emphasize that un- 
like Hod's analysis — which is an attempt to find the 
asymptotic late-time behavior in the spacetime of a spin- 
ning black hole — Poisson's analysis aims at finding the 
behavior in a spacetime in which curvature is weak ev- 
erywhere. While Poisson's analysis and results are cor- 
rect for the spacetime he studies, one should use cau- 
tion when infering from Poisson's results on the asymp- 
totic late-time behavior in a Kerr geometry: Although 
the asymptotically-far geometries are similar, the near- 
field geometries are very different. As we discuss below, 
that is a crucial element in understanding the late-time 
behavior. 

Hod's surprising predictions agree with some reported 
numerical simulations. In particular, for the case s — 0, 
£* = 0,m — Hod's formula predicts a decay rate of 
t~^, which is indeed found For the case s = 0, £* = 
4, m = 0, however. Hod's formula predicts a decay rate of 
whereas the intuitive picture predicts a decay rate of 
t^^. This case was simulated numerically by Krivan 
who found a decay rate with a non-intergal index close to 
—5.5. Like Hod, Krivan too tried to find the asymptotic 
late-time behavior in the Kerr spacetime. Some view 
this as a loose confirmation of Hod's prediction , with 
numerical accuracy of 10%, and as an invalidation of the 
intuitive picture. 

In this Rapid Communication we present results from 
independent numerical simulations for linearized pertur- 
bation fields over a Kerr background. Our simulations 
show a clear falloff rate of t^'^ for the initial data of 
s = 0, = 4, TO = 0. The quality of our results invali- 
dates Hod's prediction for the asymptotic decay rate, and 
points at difficulties with Krivan's simulations or their in- 
terpretation. In all the cases we have checked, for cither 
a scalar or a gravitational field, we find that the intuitive 
picture is correct: the late time behavior is dominated by 
the mode with the lowest value of i which can be excited. 
In particular, no spooky memory effects occur. 

We used the penetrating Teukolsky code (PTC) [l^ . 
which solves the Teukolsky equation for linearized per- 
turbations over a Kerr background in the ingoing Kerr 
coordinates {t, r, 6, (p). The Kerr metric is given by 
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cos^ 9, and M, a are the mass and 



the specific angular momentum, respectively. These co- 
ordinates are related to the Boyer-Lindquist coordinates 
{t, r, 9, ip) through ip = (p + J aA~^ dr and t = t — r + r^,, 
^ + - 2Mr and = /(r^ -|- a'^)A~'^ dr. 



where A 



const, 



Notice that t is linear in t, so that along r 
d/di=d/ dt. 

The Teukolsky equation for the function ip in the in- 
going Kerr coordinates can be obtained by implementing 
black hole perturbation theory (with a minor rescaling of 
the Kinnersley tetrad ^3)- 1^ is given by 
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Equation (0) has no singularities at the event horizon, 
and therefore is capable of evolving data across it. The 
PTC implements the numerical integration of Eq. (0) by 
decomposing it into azimuthal angular modes and evolv- 
ing each such mode using a reduced 2-1-1 dimensional 
linear partial differential equation. The results obtained 
from this code are independent of the choice of bound- 
ary conditions, because the inner boundary is typically 
placed inside the horizon, whereas the outer boundary is 
placed far enough that it has no effect on the evolution. 

The PTC has been tested in various different situa- 
tions. First, it yields the correct complex frequencies for 
the quasi-normal modes of a Kerr black hole for a wide 
range of values oia/M. Second, it has also been shown to 
yield equivalent results in the context of the close limit 
collision of two equal mass, non-spinning, non-boosted 
black holes (to ones obtained from the Zerilli formalism) 
[l3| . It is stable, and exhibits second-order convergence. 

We next set a/M = 0.9, s = 0, and T = 4, to = 0. 
The initial gaussian perturbation is taken to be a mix- 
ture of ingoing and outgoing waves, and centered about 
r = 20M with a width of 4M. As discussed above, our 
expectations are that all the even i modes are excited 
(respecting the equatorial symmetry of the initial data). 
The least £ mode which is excited is the £ = Q mode, 
so that the decay rate we expect is In contrast, 

the prediction of Hod is for a decay rate of t^^ . Fig- 
ure n shows the Teukolsky function ip for these initial 
data for = 7r/2 (the equatorial plane) for three differ- 
ent resolutions. The data clearly indicate stability and 
second-order convergence. 

A decay rate of about is already clear from Fig. ^ 
Evaluating the decay rate from the slope of the field is 
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FIG. 1: The Teukolsky function i/; as a function of t for three 
different grid resolutions, for s = 0, a/M = 0.9, £* — 4, and 
m = 0. Dotted line: 8000 steps in r and 40 steps in 8. Dashed 
line: 10000 steps in r and 50 steps in 0. Solid line: 12000 steps 
in r and 60 steps in 0. The time step is always taken to be 
half the step in r. The data are shown along r = 20M on the 
equatorial plane. 



very inaccurate: The slope then depends on the inter- 
val one chooses, and also on the presence of subdomi- 
nant modes. The first difficulty can be handled by con- 
sidering the local power index n which we define as 
n = dfip. The second difficulty can be handled by 

extrapolating n to timelike infinity. Figure |2J'\. shows n 
as a function of M/t. Timelike infinity is at zero, and 
both the regime where the field is dominated by the 
quasi-normal ringing and the regime where the field is 
dominated by the power-law tails are shown. The local 
power index n = 2.9846 at < = 1500M. Figure shows 
the behavior of 3 — n as a function of M/t. Clearly, n 
gets closer with time to the expected value of 3. In fact, 
extrapolating n to t oo using Richardson's deferred 
approach to the limit, we find the asymptotic value of n 
to be rioo = 3.0003 ± 0.0011. Our results suggest that 
the late-time field is dominated by the £ = mode. We 
checked this by plotting V' as a function of 6 in Fig. O for 
different values of i. We indeed find that ip quickly loses 
any dependence on 9, such that at late times it is indeed 
described by the £ — mode. Any dependence of -0 on 
is smaller than 3 parts in 10^ at t = lOGOAf. 

Next, we present results for the behavior of fields with 
higher spins. We set the parameters to s = 2, a/M — 0.3, 
and initial I* — 6, m = 0. The pulse is again centered 
about r = 20M with a width of AM. The prediction 
of Hod's formula for this case is a decay rate of t^^. In 
this case our expectations are that the least £ mode to be 
excited is the £ — 2 mode. Consequently, we expect the 
decay rate to be t^'^ . This is indeed confirmed in Fig.0J\, 
which shows the local power index n as a function of t/M, 
and in Fig.0}3 which displays 7 — n as a function of M/t. 
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FIG. 2: The local power index for the same data as in Fig.Q 
Upper panel (A): n as a function of M/t. Lower panel (B): 
logjQ |3 — n| as a function of M/t. The data are shown for an 
equatorial curve at r = 20M. 
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FIG. 3: The normalized Teukolsky function 'tp/%l}{6 — 0) as 
a function of at r = 25M, for the same data as in Fig. Q 
for various values of t. Upper panel: At t = 150M (dotted 
line), 200M (dash-dotted), 250M (dashed), and 300M (solid 
line). Lower panel: At f = 400M (□), 500Af (*), 600M (o), 
700 Af (-I-), 800A/ (dotted line), 900 Af (dash-dotted), 1000 Af 
(dashed), and 1500A/ (solid line). 



At t = 1500M, we find that n = 6.8646. Extrapolating 
n to timelike infinity, we find that rioo = 7.01 ± 0.03, in 
agreement with our expectations. 

Our results clearly show that starting with a pure mode 
£*,m, the late-time decay rate is dominated by the least 
mode £niin which is consistent with the equatorial sym- 
metry of the initial data and is equal to or greater than 
the least radiative mode £o — max(|s|,|TO|). The late- 
time decay rate is given by i-(2^min+3)^ q^j. conclusions 
are in sharp disagreement with the recent predictions by 
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FIG. 4: The local power index for s = 2, a/M = 0.3, T = 6, 
and m = 0. Upper panel (A): n as a function of t/M . Lower 
panel (B): log]^Q(7 — n) as a function of M /t. The data are 
taken for an equatorial curve at r = 20Af . 



Hod Hod's analysis is in the frequency domain, 

and carried to leading order in w, the angular frequency. 
That approach is very successful in the background of a 
Schwarzschild black hole, where it reproduces the known 
results . The understanding that the power-law tails 
result from scattering of the field at asymptotically large 
distances implies that it is only the small lo which are 
responsible for the tails. That is indeed the case with a 
Schwarzschild black hole. We conjecture that it would 
also be the case for a Kerr black hole, if there were no 
excitations of dominating modes which are not present 
in the initial data. For example, in the case of a scalar 
field with = m = 0, the dominating mode is already 
present in the initial data. Considering only the small 
uj contributions indeed produces a result in agreement 
with numerical simulations. When the dominating mode 
is not present in the initial data, however, it needs first to 
be excited. If it is excited (with any nonzero amplitude), 
the small w approximation may produce the correct re- 
sult for the decay rate. However, mode-excitation is an 
effect which is nonlinear in the gravitational potentials, 
and is strongest in the near zone. This suggests to us 
that a leading order (in lu) analysis will not, in general, 
get all the excited modes right. It might be the case that 
higher orders in w are necessary in order to get all the 
modes which are excited. Our numerical results indeed 
show, that when the least mode which can be excited 
^min is "far" from the initial f , that technique does not 
produce the former: For example, for initial — A and 
m — Q, the leading order in uj analysis was able to get 
the i — 1 mode excited (as is manifested by Hod's decay 
rate of t^^), but not the (. — taun — mode (which im- 
plies a decay rate of t~^). We suggest, that although a 
frequency-domain analysis is capable of getting the de- 



cay rate right, it should include an expansion to higher 
orders in to. Such an expansion would be a formidable 
endeavor. In a similar way, by taking spacetime to be 
weakly curved everywhere, Poisson tacitly assumed that 
it is just the far-zone part of the field which is important. 
(In Poisson's case, we emphasize, this assumption is well 
justified, because in the spacetime studied by Poisson 
spacetime is nowhere strongly curved. Incidentally, Pois- 
son suggests a selection-rule mechanism in the spacetime 
he studied, which is related to the remarkable vanishing 
of terms in the initial data in the transformation from 
spheroidal to spherical coordinates. The mechanism sug- 
gested by Poisson demonstrates how indeed Hod's results 
could be correct in that context. However, no such mech- 
anism is offered for a Kerr spacetime.) That assumption 
is equivalent to taking the large-r approximation, or the 
small-LJ approximation. Consequently, Poisson and Hod 
make, in fact, the same kind of approximation, such that 
it is not surprising that they obtain the same results. 
We emphasize, that Poisson acknowledges that effects 
which are nonlinear in the gravitational potentials may 
produce modes with i values which are smaller than those 
obtained by him. Poisson then remarks, that no such ef- 
fects have been reported on in the literature. Evidence 
for such an effect is precisely what we find here. Although 
the late-time expansion method does not seem to suf- 
fer from similar weaknesses, it is hard to apply for the 
problem of interest. Starting with an initial V which is 
"far" from the least mode ^min to be excited, the method 
of Ref. requires many iterations in order to find the 
excited mode i'min- Specifically, three iterations are re- 
quired in order to find the t — Q mode starting with 
£* ~ 4:,m ~ 0. Carrying this iterative scheme in practice 
seems like a daunting task. We would like to repeat, that 
while Hod's method fails to obtain the correct asymptotic 
decay rate, it may still be useful in determining an inter- 
mediate behavior for carefully chosen parameters. 

Lastly, our results are in disagreement also with the 
numerical results of Krivan [lH , who reported on a frac- 
tional power-law index which is about —5.5 for the case of 
initial s = 0, £* = A and m = 0. While we cannot point 
with certainty to the reason why Krivan's simulations 
produce a result for the asymptotic late-time behavior 
which is at odds with ours, we would like to mention 
some of the factors which may be responsible: Krivan 
takes the black hole to spin exceedingly fast. In fact, Kri- 
van takes a/M = 0.9999. The high spin of the black hole 
may act in two ways: First, it slows down the decay rate 
of the quasi-normal ringing, such that longer integration 
times are required in order to obtain the tails. Second, 
the numerical solution of the Teukolsky equation is more 
sensitive and harder when the spin is very high. Another 
factor is related to the location and the direction of Kri- 
van's initial perturbation. Krivan takes the perturbation 
to be centered around r»/Af = 100, and to have a very 
large width (of lOOM). Also, the perturbation is purely 
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outgoing on the initial slice. We thus conjecture that the 
dominating £ = mode is excited only with a very low 
amplitude, because most of the perturbation field does 
not probe the strong- field region. This, in addition to 
the great distance and width of the initial perturbation, 
may combine into late-time tails whose asymptotic be- 
havior becomes evident only at very late times, to which 
Krivan's simulations have not arrived. 

The picture which arises for linearized perturbations 
in the background of a spinning black hole is simpler 
than that which is implied by. However, we expect the 
picture to be even simpler than that for fully nonlinear 
perturbations: When the initial perturbation is not axi- 
ally symmetric, the evolving spacetime will not be axially 
symmetric either. Consequently, the m value of the field 
will not be conserved, and different values of m will also 
be excited, preserving only the equatorial symmetry of 
the initial data. We therefore expect a fully nonlinear 
evolution to yield results which arc simpler than those 
obtained from a linearized analysis: Because m is no 
longer fixed, the restriction of £o is no longer so strict: 
£q = \s\, and the dominating mode is simply the least 
£ mode which is consistent with the equatorial symme- 
try which is equal to or greater than Iq- We thus expect 
generic tails to always have a decay rate of f^f^l^l+s) . The 
more complicated results of this Rapid Communication 
then are an artifact of the linearization: the full theory 
is simpler. 
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